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Abstract 

Low temperatures may cause severe growth inhibition and mortality in fish. In order to understand the mechanism of cold 
tolerance, a transgenic zebrafish Jg (smydV.mSck) model was established to study the effect of energy homeostasis during 
cold stress. The muscle-specific promoter Smydl was used to express the carp muscle form III of creatine kinase {M3-CK), 
which maintained enzymatic activity at a relatively low temperature, in zebrafish skeletal muscle. In situ hybridization 
showed that M3-CK was expressed strongly in the skeletal muscle. When exposed to 13°C, Tg {smydl :m3ck) fish maintained 
their swimming behavior, while the wild-type could not. Energy measurements showed that the concentration of ATP 
increased in Tg {smydl :m3ck) versus wild-type fish at 28°C. After 2 h at 13°C, ATP concentrations were 2.16-fold higher in Tg 
{smydl :m3ck) than in wild-type (P<0.05). At 13°C, the ATP concentration in Tg {smydl :m3ck) fish and wild-type fish was 
63.3% and 20.0%, respectively, of that in wild-type fish at 28°C. Microarray analysis revealed differential expression of 1249 
transcripts in Tg {smydl:m3ck) versus wild-type fish under cold stress. Biological processes that were significantly 
overrepresented in this group included circadian rhythm, energy metabolism, lipid transport, and metabolism. These results 
are clues to understanding the mechanisms underlying temperature acclimation in fish. 
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Introduction 

Water temperature is considered a major environmental factor 
influencing the development, growth, reproduction, and behavior 
of aquatic ectotherms [1,2]. Fish, one kind of aquatic ectotherm, 
may face a large range of daily or seasonal temperature 
fluctuations [3]. One way to cope with these challenges is to 
migrate to areas with more suitable temperatures [4] . Another way 
is to evolve the ability to acclimate to environmental temperature 
fluctuations. Many species have evolved a variety of mechanisms 
to alter cellular biochemistry, such as by producing temperature- 
specific isoenzymes [5] , changing the membrane lipid content and 
the degree of fatty acid unsaturation [6,7], recruiting different 
muscle fiber types [8], inducing molecular chaperones [9], and 
generating antifreeze proteins during long-term evolutionary 
adaptation to a freezing environment [10]. These provide clues 
to the mechanism by which fish cope with temperature changes 
and enhance their cold tolerance. 

These "thermal compensations" may rely on energy metabo- 
lism to carry out every kind of physical activity. Muscle occupies 
more than 75% of the whole body. Muscle contraction is the main 
reaction by which chemical energy is converted to mechanical 
work and heat, the latter being essential for maintaining body 
temperature [11,12]. Muscle contraction requires large amounts of 
ATP, most of which is derived from mitochondrial electron 
transport and oxidative phosphorylation of substrates from fatty 



acid or glucose metabolism. When high-energy supplementation is 
needed, ATP can be complemented by phosphocreatine generated 
by creatine kinase. Creatine kinase (CK, EC 2.7.3.2) belongs to the 
phosphagen kinase family and functions in chemical energy 
storage and quick release; it is thus important for energy 
metabolism in muscle. CK removes a phosphate from ATP to 
produce ADP and phosphocreatine; in the reverse reaction, it adds 
a phosphate to ADP to produce creatine and ATP [13,14]. In 
common carp {Cyprinus carpio), three sub-isoforms of muscle 
creatine kinase have been designated Ml-CK, M2-CK, and M3- 
CK [15]. The formation of distinct homo- or heterodimers of M- 
CKs was altered in response to temperature acclimatization. At 
high temperature (30°C), MlMl-CK and M2M2-CK homodi- 
mers and M1M3-CK heterodimers predominate, whereas M3M3- 
CK homodimerizes at low temperature (10°C). M3-CK remained 
stable; it maintained its high enzyme activity and it was thought to 
be the dominant isoform mediating carp acclimation to lower 
environmental temperature [16]. 

The zebrafish [Danio rerio) has become one of the most 
important model organisms in many fields of research including 
genetics, neuroscience, development, physiology, toxicology, and 
biomedicine [17,18]. Zebrafish are a tropical teleost species; the 
extreme temperatures recorded in known zebrafish habitats were 
16.5°C and 38.6°C [19]. A number of other attributes have 
contributed to the scientific importance to this species, such as its 
small size, rapid development, and well-known genetic resources. 
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Zebrafish were therefore considered suitable for use as a model for 
understanding the mechanisms of acclimation under low temper- 
atures. 

The purpose of this study was to investigate the mechanism of 
cold tolerance in a transgenic zebrafish model expressing carp M3- 
CK driven by the zebrafish muscle-specific smydl promoter. Cold 
tolerance ability in the transgenic zebrafish Tg [smydl :m3ck) was 
evaluated by measuring swimming performance, balance, orien- 
tation at low temperatures, and energy production (ATP 
concentration). Microarray-based expression profiling was used 
to compare the whole transcriptome of Tg- [smydl :m3ck) and wild- 
type fish after acute low-temperature treatment. A group of genes 
related to circadian rhythm, energy metabolism, lipid transport, 
and other metabolic pathways were significantly up-regulated 
under cold-stress in the transgenic fish. These data provided 
functional genomic and physiological evidence for cold acclima- 
tion in fish. 

Materials and Methods 

Ethics statement 

Experiments involving zebrafish were conducted according to 
the regulations of local and central government, and were 
approved by the Institutional Animal Care and Use Committee 
of the Institute of Oceanology, Chinese Academy of Sciences. 

Zebrafish maintenance 

The adult captive-bred strain of zebrafish (wild-type) was 
originally obtained from an Aqua Hobby in Nanshan market 
(Qingdao, China) and acclimated to laboratory conditions. Wild- 
type and transgenic fish were reared at 28°C±0.5°C and under a 
photoperiod regime of 14/10 h light/ dark. Water quality was 
maintained by circulation with a filtration system. 

Transgenic plasmid construction 

The design of the transgenic cassette is illustrated in Figure lA. 
The smydl -gfp plasmid [20] containing the 5.3-kb smydl 
promoter fragment was used as the backbone of the transgenic 
plasmid; the gfp cassette was released by BamYLl digestion. The 
1150-bp carp mJck gene was amplified from the cDNA of 
common carp by PCR with m3ck-F and m3ck-R primers, 
according to the archived sequence (GenBank: AF055290.1); the 
15-base homologous sequences on each side of the smydl -gfp 
Bamlil site were added to the 5' end of each primer (Table 1). 
The m3ck cassette was recombined into the BamHl site of smydl - 
gfp (CloneEZ PCR Cloning Kit, GenScript, USA) to obtain the 
smydl -m3ck construct. The iTol2 cassette [21] was PCR-amplified 
from piTol2-Kan [22] with primers iTol2-F and iTol2-R, in which 
a Sail site was incorporated at the 5' end (Table 1). The PCR 
product was cloned into the Sail site of smydl -m3ck to generate 
the transgenic plasmid iTol2-smydl-m3ck. PCR cycling conditions 
for m3ck and iTol2 were as follows: 5 min initial denaturation at 
94°C; 30 cycles of denaturation at 94°C for 30 s, annealing at 
55°C for 30 s, and extension at 72°C for 1.5 min; and a final 10- 
min extension at 72°C. 

Generation of transgenic zebrafish 

Tol2 transposase RNA was synthesized in vitro using pSP64 
poly(A)/To/2 transposase plasmid [23]. The plasmid was linear- 
ized with Xbal, and capped mRNA was synthesized using 
mMESSAGE mMachine SP6 Kit (Ambion, USA). iTol2-smydl- 
m3ck plasmid DNA (50 ng/|iL) and Tol2 transposase mRNA 
(20 ng/|iL) were co-injected into one-cell stage zebrafish embryos 
(wild-type) with PLI-100 Basic Pico-Injector (Harvard, USA). 



fTo/2 cassette smyd"/ promoter m3ck 




B 

+ - + wt G P B 



m3ck 
HH 




Figure 1. The transgenic cassette and identification of 
transgenic zebrafish. A. Design of the transgenic cassette. mSck 
was driven by tlie sl<eletal muscle-specific promoter smydl and tlie 
/To/2 cassette was introduced for germline transgenes. B. Identification 
of transgenic zebrafish by PCR with m3ck primers. HH-F/R primers 
derived from tiggy winkle hedgehog exon 3 were used as control 
primers. non-transgenic fish; +, transgenic fish; wt, embryo of wild- 
type fish; G, genomic DNA of wild-type zebrafish; P, the transgenic 
plasmid iTol2-smyd1-m3ck as positive control; B, blank control. 
doi:1 0.1 371/journal.pone.01 02492.g001 

Injected fish were cultured to sexual maturity and outcrossed to 
identify germline transgenic carriers. 

Identification of transgenic zebrafish 

To identify FO transgenic zebrafish, each putative founder was 
reared to sexual maturity and outcrossed with wild-type fish to 
obtain the corresponding Fl progenies. For each founder, 100 Fl 
embryos were collected at 72 h post-fertilization (hpf) and 
anesthetized with 100 mg/L MS222 (Tricaine methane sulfonate). 
The embryos were submerged in 500 |iL genomic DNA extraction 
buffer (10 mM Tris pH 8, 2 mM EDTA, 0.2% Triton X-100, 
200 |lg/mL Proteinase K) and ground thoroughly. The embryos 
were incubated at 50°C and vortexed occasionally until completely 
dissolved. After boiling in a water bath for 10 min, samples were 
centrifuged at 12,000 rpm for 5 min; the supernatant was stored at 
— 20°C. PCR was performed in a total volume of 25 |LiL with 
1.5 |iL supernatant. The m3ck-F/R primers were used to detect 
transgenic embryos. The HH-F/R primers derived from the tiggy 
winkle hedgehog exon 3 sequence were used for the control [24] 
(Table 1) with the cycling conditions described above. The 
transgenic zebrafish were named Tg [smydl :m3ck). 

The Fl generation was produced by outcrossing FO with wild- 
type; the F2 generation was produced by an inter-cross of Fl 
progeny. To identify transgenic zebrafish of the Fl and F2 
generation, fish were cultured to age two months. Half of the tail 
fin was cut and sampled as above using 40 |lIL of genomic DNA 
extraction buffer. PCR was performed as described above. 

Whole-mount in situ hybridization 

To detect muscle-specific m3ck expression, whole-mount in situ 
hybridization was performed on the corresponding Fl transgenic 
fish. The FO Tg [smydl :m3ck) was crossed with wild-type to 
produce an Fl generation. The Fl embryos were raised to 72 hpf 
and anesthetized with 100 mg/L MS222, fixed in 4% parafor- 
maldehyde/PBS overnight at 4°C, and bleached with bleaching 
buffer [according to the Zebrafish Book, 5th Edition, Westerfield 
(2007)]. The embryos were washed in PBST, transferred to 
methanol, and stored at — 20°C. Whole-mount in situ hybridiza- 
tion with a digoxigenin (DIG)-labeled antisense RNA probe was 
performed as described [25]. The carp m3ck PCR product was 
cloned into the pBluescript II SK (+) Sail site to generate the 
m3ck-BSSK fusion construct. The m3ck-BSSK plasmid was 
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Table 1. PCR primers for constructing and identifying transgenic zebrafish. 



Primer 


Sequence 


m3ck-F 


5 ' -AAGAATTCCTCG ACGCATG ACTAAG AACTG-3 ' 


m3ck-R 


5 ' - GCTCTAG AACTAGTGGCCATTTTACTTCTG-3 ' 


iTol2-F 


5 '-GTCGACCCCTGCTCG AGCCGGGCCCAAGTG-3 ' 


iTol2-R 


5 '-GTCGACATTATGATCCTCTAGAT CAGATCT-3 ' 


HH-F 


5'- GGACGGTGACACTTGGTGATG-3 


HH-R 


5'- CGAGTGGATGGAAAGAGTCTC-3 


Tol2L-out 


5 '-CCCTGCTCG AGCCGGGCCCAAGTG-3 ' 


Tol2R-out 


5 '-ATTATGATC-CTCTAGATCAGATCT-3 ' 



Homologous sequence was highlighted in bold. 
doi:1 0.1 371/journal.pone.01 02492.t001 



linearized with EcoKl and used to synthesize the riboprobe by 
using a DIG RNA Labeling Kit (Roche, Germany) by in vitro 
transcription with T3 RNA polymerase. Immunoreactive signal 
was detected by color microscopy. 

Gene insertion site analysis 

Inverse PCR was performed to analyze the insertion site in the 
zebrafish genome. Genomic DNA from transgenic zebrafish was 
extracted and 3 |ig was digested with Alul or Haelll for 5 h. After 
self-ligation with T4 Ligase (Promega, USA) for 24 h at 15°C, 
circular genomic DNA was obtained and used as the PCR 
template in a 50-|iL reaction containing 200 ng DNA, 0.4 |iM 
Tol2L-out and Tol2R-out primers [21] (Table 1), and 25 |iL of 
2x Taq Polymerase Mix (CWBIO, China). Cycling conditions 
were as follows: 5 min initial denaturation at 94° C; 35 cycles of 
denaturation at 94°C for 30 s, annealing at 52 °C for 30 s and 
extension at 72°C for 2 min; and a 10 min final extension at 72°C. 
The PCR product was sequenced and compared to the zebrafish 
genome database to identify the insertion sites. 

Cold stress 

The F2 generation was raised to the age of 3 months 
(2. 17 ±0.02 cm, 0.42 ±0.01 g) before use in the cold acclamation 
experiment. Consistent with previous reports [26], 13°C was 
selected to represent cold stress. Treatment time was determined 
according to the (1) experimental feasibility; and (2) time required 
to discern a differential response between groups. Cold stress was 
applied for 2 h according to Long et al. [27] and our own 
preliminary experiments. Before cold treatment, the fish were 
transferred to tanks with small openings in the bottom and allowed 
to stand for one day. They were gently moved to the cold bath for 
the experiment. The groups (n = 30 per group) were designated 
m3ck-13°C, m3ck-28°C, wt-13°C, and wt-28°C. After tempera- 
ture stress for 2 h, the fish were anesthetized with 100 mg/L 
MS222. To determine the ATP concentrations, fish were dissected 
on ice and the muscle was sampled. For microarray or real-time 
quantitative PCR analysis, the whole fish was transferred 
immediately into liquid nitrogen. 

Deternnination of ATP concentration 

After cold stress, five fish from each group were sampled. 
Intracellular ATP levels were determined using a bioluminescence 
ATP assay kit (Beyotime, China). Approximately 30-40 mg 
muscle samples were removed from each fish, homogenized in 
lysis buffer at a final concentration of 10 mg/100 |iL, and 
centrifuged at 1 2,000 xg to collect the supernatant. ATP detection 



working solution (100 |iL) was added to each well of a 96-well 
white culture plate and incubated for 3-5 min at room temper- 
ature. Then 80 |iL of lysis muscle solution and 100 |iL of 
luciferin-luciferase reagent were mixed for 3 s before lumines- 
cence was measured over 10 s on an Infinite F200 PRO 
multimode reader (TEC AN, USA). Each lysate was examined in 
three parallel wells. The 96-well plates also contained wells with 
serial dilutions of an ATP standard to generate a standard curve. 
Results were expressed as relative ATP levels. 

Microarray analysis 

Gene expression in triplicate samples of m3ck-13°C, m3ck- 
28°C, wt-13°C, and wt-28°C was assessed by Agilent Zebrafish 
Oligo Microarray 4x44 K (V3) (design ID: 026437) at Shang- 
haiBio Corporation (SBC). Twelve microarray experiments were 
performed, each with three fish. When encountering cold stress, 
fish perform as an integrated system and cold tolerance is judged 
from the whole-body response. So, RNA samples were extracted 
from whole zebrafish with Trizol Reagent (Invitrogen, CA, USA). 
RNA quality as expressed by the RNA integrity number (RIN) was 
determined on an Agilent Bioanalyzer 2100 (Agilent Technologies, 
Santa Clara, CA, USA) and then purified with an RNeasy mini kit 
(QIAGEN GmBH, Germany) with the RNase-Free DNase Set 
(QIAGEN). Total RNA was amplified and Cy3-CTP labeled with 
a Low Input Quick Amp Labeling Kit, One-Color (Agilent). 
Labeled cRNA was purified with the RNeasy mini kit (QIAGEN). 
Slides were hybridized with 1.65 |ig Cy3-labeled cRNA, pro- 
cessed, and scanned on an Agilent Microarray Scanner (Agilent) 
with default settings. Data were extracted with Feature Extraction 
software 10.7 (Agilent). Raw data were normalized using the 
Quantile algorithm in Gene Spring Software 11.0 (Agilent). The 
raw and normalized datas were deposited in the NCBI's Gene 
Expression Omnibus (GEO) database (Accession NO. GSE58988). 
Differentially expressed gene identification (FC >2, p<0.05) and 
pathway enrichment analysis were performed by the SAS online 
system of ShanghaiBio Corporation. GO enrichment analysis was 
performed with the SAS online system and the GOEAST web- 
based software (http://omicslab.genetics.ac.cn/GOEAST/index. 
php). The Venn diagram was drawn with R package "VennDia- 
gram". 

Real-tinne quantitative PCR 

To validate the microarray data, expression of per lb, nrldl, 
per2, dbp2, opnlmw2, tmx3, and ascllh was determined by qPCR 
analysis. These genes were selected for their microarray-derived 
expression patterns and based on previous reports on cold stress 
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response. Samples were analyzed in two batches. One batch was 
the same as that used in the microarray analysis. The second batch 
included samples from the timed cold treatments (five fish were 
randomly sampled from each group). 

First-strand cDNA was synthesized from 2 jig total RNA using 
random hexamer primers. PGR primers were designed using 
Primer Premier 5.0 software (Premier, Canada). qPCR was 
performed in a Mastercycler ep realplex (Eppendorf). The 20-|iL 
reactions contained 2 x SYBR Green Real-Time PGR Master mix 
(TAKARA), 2 pmol each primer, and 10 x diluted cDNA 
templates. All reactions were performed in triplicate with the 
following cycling conditions: 40 cycles of 5 s at 95°G, 30 s at 56°G 
and 20 s at 72°G, followed by melt curve analysis. Reaction 
specificity was confirmed by the observation of a single melt peak 
at the expected Tm. The sequences, Tm, and amplification 
efficiency, accession number, gene name, and amplicon lengths 
are listed in Table SI. 

Statistical analysis 

Swimming status and ATP concentrations were analyzed by 
one-way ANOVA followed by ^test. qPGR results were analyzed 
by two-way ANOVA followed by Bonferroni post-tests. Significant 
difference was defined as P<0.05. Data represent mean ± S.E.M. 

Results 

A transgenic line expressing IVIBCK in skeletal muscle 

One out of nine putative founders was found to carry the 
smydl:m3ck transgene by germ-line transmission and selected to 
generate the transgenic line (Figure IB). In situ hybridization 
confirmed that carp m3ck was specifically expressed in the skeletal 
muscle (Figure 2). The carrying rates of the Fl and F2 generation 
were 16.22% (6/37) and 34.9% (22/63), respectively. The 
transgenic cassette was proved to be inserted into a non-coding 
area by inverse PGR. 



A 



B 




Figure 2. Expression of m3ck in transgenic fish by in situ 
hybridization. Whole-mount in situ hybridization was used to detect 
mic/c expression in wild-type (A) and Tg [smydlimScK) fish (B) at 72 hpf. 
The signal was evaluated by microscopy for color detection of 
immunoreactive signals. Scale bars = |am. 
doi:l 0.1 371 /journal. pone.0102492.g002 



SwimiTiing behavior and energy 

During exposure to 13°G water temperature, the wild-type fish 
instantly lost their swimming ability. They did not move; they were 
imbalanced, inverted, disoriented, and did not respond to external 
stimulation such as gentle knocking on the side of the container. In 
contrast, Tg [smydl:m3ck) fish did not move or were inverted for 
the first 5-10 min, after which normal swimming was restored 
(Video SI and Video S2). The percentages of fish exhibiting 
smooth swimming, occasional swing, and response to external 
stimulation were 23.3% (7/30), 73.3% (22/30), and 70.0% (21/ 
30), respectively, in Tg {smydl:m3ck) fish and 0, 36.7% (11/30), 
and 20.0% (6/30) in wild-type (P<0.05). The percentages of 
transgenic fish exhibiting no movement [3.3% (1/30)] and 
inversion [6.7% (2/30)] were lower than in wild-type fish 
[63.3% (19/30) and 23.3% (7/30); P<0.05]. The number of 
deaths did not significantly differ between groups (Figure 3). Thus, 
we concluded that our transgenic fish performed better than wild- 
type under cold stress. 

ATP concentrations were 10.04 |iM at 28°G and 2.01 |iM at 
13°G in wild-type fish and 12.66 |iM at 28°G and 6.36 |iM at 
13°G in Tg {smydl:m3ck) fish (Figure 4). Thus, the concentration 
of ATP tended to increase in Tg {smydl:m3ck) versus wild-type fish 
at 28°G. After 2 h at 13°G, the ATP concentration in Tg 
[smydl:m3ck) fish was 2.16-fold greater than in wild-type fish {P< 
0.05). If we define the wild-type fish at 28°G as normal, the ATP 
concentration at 13°G was 63.3% in Tg {smydl:m3ck) fish, but 
only 20.0% in wild-type fish. 

Microarray analysis 

Agilent microarray analysis was used to detect whole-body gene 
expression in transgenic and wild-type fish. Principle Gomponent 
Analysis (PGA) was performed to verify consistency within groups 
and differences between groups. There were 507 up-regulated 
transcripts and 742 down-regulated transcripts in the m3ck-13°G 
vs. m3ck-28°G comparison. In the wt-13°G vs. wt-28°G compar- 
ison, there were 165 up-regulated transcripts and 105 down- 
regulated transcripts (Figure 5). Under cold stress, there were 210 
up-regulated and 252 down-regulated transcripts in Tg 
[smydT.m3ck) vs. wild-type fish and at 28°G, there were 1095 
up-regulated transcripts and 466 down-regulated transcripts in 
m3ck-28°G vs. wt-28°G fish. 

Validation by quantitative real-time PCR 

To validate the microarray results, relative mRNA levels for 
seven genes were measured by qPGR (Figure SI). The results 
showed excellent agreement with the microarray results and 
confirmed the reliability of the single-channel microarray analysis 
method. 

Overlapping expressed genes under cold stress 

After 2 h under cold stress, more genes were differentially 
expressed genes in Tg {smydl:m3ck) than in wild- type fish 
(Figure 6). Thirty up-regulated transcripts were shared by m3ck- 
13°G vs. m3ck-28°G and wt-13°G vs. wt-28°G (Figure 6A), 
corresponding to cry lb (cryptochrome lb), cry 3 (cryptochrome 5), 
cry 3 (cryptochrome 3), nos2a (nitric oxide synthase 2a, inducible), 
illh (interleukin 1, het3), per lb [period homolog lb {Drosophila)], 
per2 [period homolog 2 {Drosophila)] , egrl (early growth response 
1), nrldl (nuclear receptor subfamily 1, group d, member 1), dbp2 
(D-box binding protein 2), and slc35f3 (solute carrier family 35, 
member F3). Twenty-eight genes were down-regulated (Figure 6B), 
including opnlmwl [(opsin 1 (cone pigments), medium-wave- 
sensitive I], opnllw2 [opsin 1 (cone pigments), long-wave-sensitive 
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Figure 3. Status of transgenic and non-transgenic fish at 13°C The percentages smooth swimming, occasional swing, response to external 
stimulation, no movement, inversion, and dead fish were recorded. Data were analyzed by one-way ANOVA and t-test to compare data between fish 
lines. Data represent mean ± S.E.M. (n = 30). '""'indicates P<0.05. 
doi:1 0.1 371 /journal. pone.01 02492.g003 



2], tmx3 (thioredoxin-related transmembrane protein 3), grklb (G 
protein-coupled receptor kinase 1 b), nmnatl (nicotinamide 
nucleotide adenylyltransferase 1), trpv6 (transient receptor poten- 
tial cation channel, subfamily V, member 6), gucale (guanylate 
cyclase activator le), taarl4h (trace amine associated receptor 
14 h), itpr3 (inositol 1 ,4,5-triphosphate receptor, type 3), abcg4b 
[ATP-binding cassette, sub-family G (WHITE), member 4b], and 
ascllb [achaete-scute complex-like lb [Drosophila)]. However, 
only one transcript was down-regulated in m3ck-13°C vs. m3ck- 
28°C but up-regulated in wt-13°C vs. wt-28°C: retinol dehydro- 
genase 13 -like (Figure 6G). No transcript was found to be up- 



regulated in m3ck-13°C vs. m3ck-28°C and down-regulated in wt- 
13°C vs. wt-28°C (Figure 6D). 

Gene Ontology (GO) enrichment analysis 

To identify transcriptional events that occur during cold 
acclimatization, up- and down-regulated genes in each group 
were subjected to GO enrichment analysis (Table S2). In wild-type 
fish, cold stress induced biological processes such as regulation of 
cellular processes, regulation of DNA-dependent transcription, 
DNA metabolic process, response to stress, and protein-DNA 
complex subunit organization (Figure 7 A and Figure S2). Up- 

wild-type 
Tg (smyd1: m3ck) 



0. 

I- 
< 

o 

c 
o 

c 
o 
o 
c 
o 
o 



20- 
15- 
10- 
5- 
0- 



28 X 13 X 

Temperature 



Figure 4. ATP concentration in transgenic and non-transgenic fish at 28°C or 13°C The concentration of ATP in Tg {smydhmSck) and wild- 
type fish was measured at the indicated times. Data were analyzed by one-way ANOVA and t-test. Data represent mean ± S.E.M. {n = 5). "^"indicates 
P<0.05. 

doi:1 0.1 371 /journal.pone.01 02492.g004 
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Figure 5. Relative abundance of up- and down-regulated 
transcripts in the experimental groups. The numbers of differen- 
tially expressed genes are shown. Transcripts with fold change >2, 
P value <0.05, were considered to be differentially expressed. 
doi:1 0.1 371 /journal. pone.01 02492.g005 

regulated molecular functions included nucleic acid binding 
transcription factor activity, signal transducer activity, lipid 
binding, and thyroid hormone receptor activity (Figure 7B). In 
transgenic fish, cold-induced biological processes included lipid 
transport, DNA metabolic process, responses to endogenous 
stimulus, membrane organization, and clustering of voltage-gated 
sodium channels (Figure 7C). Among these processes, lipid 
transport was the most enriched term (Figure S3). Up-regulated 
molecular functions included lipid transporter activity, substrate- 
specific transporter activity, and acyl-CoA desaturase activity 
(Figure 7D). 

Biological processes such as organic acid metabolic process, 
carboxylic acid metabolic process, and oxoacid metabolic process, 
were enriched among m3ck-13°C vs. wt-13°C cold stress-induced 
genes (Figure 8A). At normal temperature, however, overrepre- 
sented GO terms among m3ck-28°C vs. wt-28°G up-regulated 



genes were peptidyl-histidine phosphorylation, oligopeptide trans- 
port, and meiosis I (Figure 8B). 

Biological processes inhibited by cold stress in wild-type fish 
included the G-protein coupled receptor signaling pathway, 
sensory perception of light stimulus, protein-chromophore linkage, 
and dicarboxylic acid metabolic process (Figure S4). Molecular 
functions included G-protein coupled receptor activity, iron ion 
binding, photoreceptor activity, malate dehydrogenase activity, 
and biotin binding (Figure 9A). For transgenic fish, cold stress- 
inhibited biological processes included the G-protein coupled 
receptor signaling pathway, sensory perception of light stimulus, 
protein-chromophore linkage, and phototransduction (as in wild- 
type). Molecular functions included G-protein coupled receptor 
activity, cyclic-nucleotide phosphodiesterase activity, and guany- 
late cyclase activator activity (Figure 9B, Figure S5). Biological 
processes of the sensory perception of light stimulus and molecular 
function of calcium-sensitive guanylate cyclase activator activity 
were enriched in m3ck-13°C vs. wt-13°C cold stress-inhibited 
genes (Figure lOA). At normal temperatures, down-regulated 
biological processes in m3ck-28°C vs. wt-28°G were the establish- 
ment of localization, organic substance transport, and lipid 
transport. Molecular functions included substrate-specific trans- 
porter activity; lipid transporter activity; and transferase activity, 
transferring glycosyl groups (Figure lOB). 

Pathway enrichment in cold stress-regulated genes 

In wild-type fish, cold-induced pathways included MAPK 
signaling, circadian rhythm, cytosolic DNA-sensing, pentose and 
glucuronate interconversions, and NOD-like receptor signaling, 
toll-like receptor signaling pathway, and apoptosis. In transgenic 
fish, cold-induced pathways included the molecular clock, steroid 
biosynthesis, cyanoamino acid metabolism, taurine and hypotaur- 
ine metabolism, folate biosynthesis, and terpenoid backbone 
biosynthesis. These pathways were mainly involved in the 
regulation of biological processes such as environmental adapta- 
tion, signal transduction, membrane construction, and cell 
adhesion (Figure 1 lA, Table S3). In addition, cold stress inhibited 
primary bile acid biosynthesis, pyruvate metabolism, ether lipid 




Figure 6. Overlapping expressed transcripts under cold stress. Venn diagrams represent the number of overlapping and differentially 
expressed genes between A. Up-regulated transcripts under cold stress (2 h) of Tg {smydhmSck] and wild-type fish. B. Down-regulated transcripts 
under cold stress of Tg {smydhmSck) and wild-type fish. C. Down-regulated transcripts under cold stress of Tg {smydl :m3ck] and up-regulated 
transcripts in wild-type fish. D. Up-regulated transcripts under cold stress of Tg {smydl :m3ck) and down-regulated transcripts in wild-type fish. 
doi:1 0.1 371 /journal. pone.01 02492.g006 
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Figure 7. GO enrichment analysis of up-regulated transcripts in transgenic and wild-type fish under cold stress. A. Biological processes 
of wt-13°C vs. wt-28°C up-regulated genes. B. Molecular functions of wt-13°C vs. wt-28°C up-regulated genes. C. Biological processes of m3ck-13°C 
vs. m3ck-28°C up-regulated genes. D. Molecular functions of m3ck-13°C vs. m3ck-28°C up-regulated genes. The horizontal axis represents transcript 
count in each GO term. Up-regulated transcripts were submitted for GO enrichment analysis using the SAS online analysis system. 
doi:1 0.1 371 /journal. pone.01 02492.g007 
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Figure 8. GO enrichment analysis of up-regulated transcripts between transgenic fish and wild-type fish under cold stress or 
normal temperature. A. m3ck-28°C vs. wt-28°C up-regulated transcripts. B. m3ck-13°C vs. wt-13°C up-regulated transcripts. The horizontal axis 
represents transcript count in each GO term. Up-regulated transcripts were submitted for GO enrichment analysis using the SAS online analysis 
system. 

doi:1 0.1 371 /journal. pone.01 02492.g008 



PLOS ONE I www.plosone.org 



7 



July 2014 I Volume 9 | Issue 7 | e102492 



Cold Tolerance in Transgenic Zebrafish 



wt-13 X vs wt-28 '*C down-regulated 



Molecular function 



Biological_j)rocess 



Cellular component 



malate dehydrogenase activity 
iron ion binding 
biotin binding 
G-protein coupled receptor activity 
photoreceptor activity 
protein-chronfK)phore linkage 
sensory perception of light stimulus 
G-protein coupled receptor signaling pathway 
dicarboxylic acid metabolic process 
nonmoble primary cilium 
photoreceptor outer segment 



B 



10 



m3ck-13 °C vs m3ck-28 °C down-regulated 



Molecular function 



Biological_process 



Cellular component 



cyclic-nucleotide phosphodiesterase activity 
guanylate cyclase activator activity 
G-protein coupled receptor activity 
protein-chromophore linkage 
phototransduction 
sensory perception of light stimulus 
G-protein coupled receptor signaling 
cell projection 
nonmoiiie primary ciiium 
photoreceptor outer segment 




— 1— 

20 



— I — 
40 



— I 
60 



Figure 9. GO enrichment analysis of down-regulated transcripts of transgenic and wild-type fish under cold stress. A. wt-1 3°C vs. wt- 
28°C down-regulated transcripts. B. m3ck-13°C vs. m3ck-28°C down-regulated transcripts. The horizontal axis represents transcript count in each GO 
term. Down-regulated transcripts were submitted for GO enrichment analysis using the SAS online analysis system. 
doi:1 0.1 371 /journal. pone.01 02492.g009 
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Figure 10. GO enrichment analysis of down-regulated transcripts between transgenic and wild-type fish under cold stress or 
normal temperature. A. m3ck-28°C vs. wt-28°C down-regulated transcripts. B. m3ck-13°C vs. wt-13°C down-regulated transcripts. The horizontal 
axis represents transcript count in each GO term. Up- and down-regulated transcripts were submitted for GO enrichment analysis using the SAS 
online analysis system. 
doi:10.1371/journal.pone.0102492.g010 
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Figure 1 1 . KEGG pathway analysis of cold-regulated genes. A. Up-regulated genes. B. Down-regulated genes. Up- and down-regulated genes 
were submitted for KEGG enrichment analysis using the SAS online analysis system. Color scales represent P values of enrichment tests and gray cells 
indicate an empty value or P>0.05. 
doi:1 0.1 371 /journal. pone.01 02492.g01 1 



metabolism, and citrate cycle (TCA cycle) in wild-type fish. Some 
cold-inhibited pathways such as cytokine-cytokine receptor inter- 
action, the PPAR signaling pathway, circadian rhythm, and drug 
metabolism were overrepresented in transgenic fish (Figure IIB, 
Table S3). It is important to note that some enriched pathways or 
GO terms may not reflect actual situations in some tissues or cells 
due to our use of whole organisms for this study. 

Discussion 

In this study, we established a transgenic zebrafish line Tg 
{smydl:m3ck) with energy-producing related genes to investigate 
the mechanism of cold tolerance. Microarray analysis revealed 
numerous genes, key biological processes, and intracellular 
pathways associated with cold acclimation. These results might 
offer insights into the molecular mechanisms underlying gene 
expression during cold stress in fish. 

To survive in the wild, zebrafish must adapt to changes in water 
temperature. Daily variations of water temperature in the Ganges 
basin, a habitat of wild zebrafish, were recorded as ranging from a 
low of ±0.1 °G to a high of ±5.6°C, with seasonal variations of 
15.0±0.4°C. The lowest recorded water temperature in a 
zebrafish habitat was 16.5°C, the highest was 38.6°G [19]. Several 
factors influence thermal tolerance in fishes; however, the initial 
temperature extremity may influence tolerance of temperature 
shifts regardless of sex [28,29]. Zebrafish maintained under 
variable temperature conditions exhibit greater tolerance than 
fish maintained at constant temperature [30]. These observations 
suggest pure zebrafish lines in laboratory culture under constant 



temperature conditions (28°C) may exhibit different tolerances 
than wild zebrafish. Therefore, we used wild-type strains of 
zebrafish to study cold temperature tolerance. 

We analyzed gene expression profiles in transgenic and wild- 
type fish at 28°G to detect the impact of the transgene insertion on 
general metabolism and physiology. Pathways such as purine 
metabolism, PPAR signaling, and fructose and mannose metab- 
olism were up-regulated in transgenic versus wild-type fish. The 
up-regulation of purine metabolism may be due to physiological 
adjustments to remove excess ATP generated by the creatine 
kinase transgene. PPARs are the major regulators of lipid and fatty 
acid metabolism and regulate transport, oxidation, storage, and 
synthesis of fatty acids [31]. Moreover, PPAR activation can 
enhance sugar metabolism, which may contribute to the up- 
regulation of fructose and mannose metabolism. Pathways such as 
biosynthesis of unsaturated fatty acids, aminoacyl- tRNA, and 
other metabolic pathways were down-regulated in transgenic 
versus wild-type fish. These genes included acoxJ (acyl-coA 
oxidase 3), scd (stearoyl-CoA desaturase), and mars (methionyl- 
tRNA synthetase). AcoxJ and scd mediate fatty acid oxidation 
[32,33] and down-regulation might indicate that the overall 
energy supply in the transgenic fish was sufficient to alleviate the 
need for fatty acid metabolism. This assumption could be linked to 
up-regulation of the PPAR signaling pathway, which mediates 
fatty acid storage. Mars plays a critical role in protein biosynthesis 
by charging tRNAs with their cognate amino acids [34]. Down- 
regulation of this process in transgenic fish might also due to the 
sufficient energy supply and the reduced need for translation of 
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many metabolic enzymes. In comparison to wild-type fish, the 
transgenic fish used different metabolic strategies to maintain 
energy homeostasis. In addition, there was no difference observed 
in transgenic fish versus wild-type fish with and without cold stress 
(Figure 1 1). This would suggest that the differential pathways 
detected were not affected by the transgene manipulate. 

A common response to stimuli and the initiation of defensive 
compensatory mechanisms has been detected in transgenic and 
non-transgenic fish under cold stress. Under cold-stress, transgenic 
and non-transgenic fish down-regulate opsin 1 (cone pigments), 
opnlmwl, opnlmw4, opnlmw2, and opnllw2. These are rod cell 
photoreceptors belonging to the G-protein coupled transmem- 
brane receptor family that mediates transmission of optical signals 
to the brain [35] and thus suggest the importance of the neuronal 
signal transduction and central nervous system (GNS). It is well 
known that the CNS of both mammals and fish play important 
roles in thermoregulation. The hypothalamus is the most notable 
central thermoregulatory site to trigger the primary response to 
cold and helps to increase sympathetic output along with the spinal 
cord to peripheral tissues such as muscle and brown fat [36—38]. 

Up-regulated genes in transgenic and non-transgenic fish under 
cold stress included GO terms for DNA-dependent negative 
regulation of transcription and response to light stimulus. These 
transcripts included the cryptochrome family {cry-dash, Cry 3 and 
Cry 5) and period homolog gene family (perlb and per 2). Up- 
regulation of these genes leads to inhibition of the clock gene, the 
core of the circadian rhythm pathway, and thus down-regulation 
of the remainder of the pathway. A correlation between 
temperature and circadian rhythm has been reported in zebrafish 
[27] and mice [39]; this was considered a "temperature 
compensation" mechanism to correct the natural tendency of 
the rate of biochemical reactions to change with temperature [40]. 

Lipid transport and metabolism-related genes such as scd 
(stearoyl-GoA desaturase, delta-9-desaturase), sgms2 (sphingomy- 
elin synthase 2), fdftl (squalene synthetase), and ebp (sterol 
isomerase) were only up-regulated in transgenic fish. The 
pathways involving these genes included the biosynthesis of 
unsaturated fatty acids, sphingolipid metabolism, and steroid 
biosynthesis. Among these, scd mediates biosynthesis of unsatu- 
rated fatty acids. Delta-9 catalyzes the synthesis of monounsatu- 
rated fatty acids (18:1 n-9 and 18:2 n-6), which are required for 
triglyceride energy storage and contributes to the maintenance of 
cell membrane fluidity [41-43]. Genes such as fdftl and ebp 
participate in steroid biosynthesis and bile acid biosynthesis. The 
activity of fdftl could provide precursors to cholesterol synthesis. 
Cholesterol is an essential structural component of cell membranes 
and serves as a precursor for the synthesis of steroid hormones 
such as bile salts, which are important in fat absorption [44]. 
Sgms2 is associated with sphingolipid metabolism. Sphingolipids 
are enriched in the plasma membrane and are important 
structural lipids. They participate in the formation of membrane 
domains with cholesterol and affect the distribution and function 
of membrane proteins. Some less abundant sphingolipids are 
important modulators of cell signaling [45-47] . Thus, our results 
suggest transgenic fish could maintain their cell membrane 
construction and fluidity under lower temperatures, unlike the 
wild-type fish. Cell membrane integrity and fluidity are essential 
for material transport and cell communication. Cell membrane 
construction-related pathways such as lipid metabolism and 
primary bile acid biosynthesis were also down-regulated in wild- 
type fish under cold stress. Therefore, cell membrane fluidity may 
be important for low- temperature tolerance in fish. 



Our results demonstrated the importance of energy for fish 
during cold stress. ATP is produced by photophosphorylation, 
cellular respiration, and fermentation and is involved in many 
cellular processes, including biosynthetic reactions, motility, and 
cell division [48]. The concentration of ATP in the transgenic fish 
at low temperatures was 60% (6.36/10.04) of that of wild-type fish 
at normal temperatures, which was three times (6.36/2.01) that of 
wild-type fish at low temperatures. This confirmed that M3CK 
maintains its activity and catalyzes the formation of ATP at low 
temperatures. Some of this energy may be used for cellular 
processes such as cell membrane construction in response to low 
temperature, as indicated by the up-regulation of genes partici- 
pating in steroid biosynthesis, folate biosynthesis, and terpenoid 
backbone biosynthesis only in the transgenic fish (Figure 10). It 
should be noted that the transcripts of genes involved in oxidative 
phosphorylation, the main source of ATP, were significantly 
down-regulated in the transgenic fish. However, pathways that lie 
upstream of oxidative phosphorylation, such as pyruvate metab- 
olism and the TGA cycle, which produce substrates for oxidative 
phosphorylation, were down-regulated only in wild-type fish at low 
temperature. Therefore, oxidative phosphorylation could not 
supply enough energy in response to cold stress. It was interesting 
to note that the transcripts of genes involved in pentose and 
glucuronate inter-conversions were up-regulated in wild-type fish 
under cold stress. This is another metabolic pathway to producing 
glucose and ribose-5-phosphate (R5P), which could be used in the 
synthesis of nucleotides and nucleic acids. Therefore, glucose 
metabolism cannot be the main supplier of energy, and other 
energy-producing pathways are more important during cold stress. 
Further studies on the supply of energy should be performed. 

Transcripts of the jun and fos genes were found only in cold 
treated wild-type fish. They were also up-regulated in previous 
studies in zebrafish [49] and rats [50,51] under cold stress. Fos can 
form heterodimers with Jun, forming a transcription factor 
complex that often induces apoptosis [52,53] under stress, and is 
important for maintaining cell homeostasis. During cold stress, cell 
damage may trigger the immune response, and many genes 
involved in immune response pathways were significantly up- 
regulated in cold-treated wild-type fish, such as Toll-like receptor 
signaling [54], cytosolic DNA-sensing [55], NOD-like receptor 
signaling [56], and cytokine-cytokine receptor interactions [57] 
(Figure 10). Genes involved in the synthesis of nucleotides and 
nucleic acid-pentose and glucuronate inter-conversions were up- 
regulated in wild-type fish under cold stress (Figure 10). During the 
immune response, apoptosis functions as a defense reaction to 
remove damaged cells; however, the genes involved in cytokine- 
cytokine receptor interactions were down-regulated in transgenic 
zebrafish under cold stress. Thus, injury may cause cell death at 
low temperatures and injury might be another important concern 
for fish under cold stress. 

In conclusion, we identified multiple biological processes, 
intracellular signaling pathways, and key genes associated with 
cold tolerance by using muscle-specific transgenic zebrafish. The 
regulation of the energy metabolic process, cell membrane 
construction and fluidity, and cell death were important for 
acclimation to low temperatures. Our results also suggest energy is 
a critical factor in low-temperature adjustment. For further 
validation of our hypothesis, protein changes under cold stress 
will be a very interesting experiment that we plan to pursue in our 
future work, in addition to a study of metabolite profiles. 
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Supporting Information 

Figure SI Validation of microarray data by real-time 
quantitative PGR. qRT-1 represented qPCR with the same 
samples used in the microarray. qRT-2 represented qPCR with 
the samples cold treated at different times (n = 5). Data of qPCR 
referred to the left axis of ordinate and data of microarray analysis 
referred to the right axis of ordinate. 
(TIF) 

Figure S2 Biological processes of m3ck-13°C vs. m3ck- 
28 °C up-regulated transcripts. 

(PDF) 

Figure S3 Biological Processes of wt-13°C vs. wt-28°C 
up-regulated transcripts. 

(PDF) 

Figure S4 Biological processes of m3ck-13°C vs. m3ck- 
28 °C down-regulated transcripts. 

(PDF) 

Figure S5 Biological Processes of wt-13°C vs. wt-28°C 
down-regulated transcripts. 

(PDF) 

Table SI qRT-PCR primers used for validating micro- 
array data. 

(XLS) 

Table S2 GO enrichment analysis of cold-regulated 
genes. 

(XLS) 



Table S3 Pathways enrichment analysis of cold-regu- 
lated genes. 

(XLS) 

Video SI Swimming behavior of wild-type zebrafish at 
13°C. The timeline represents: 00:00-00:27: After the fish moved 
into cold water bath and the beginning of cold treatment; 00:28— 
01:02: after cold treated for 1 h; 01:03-01:36: after cold treated for 
2 h. 
(MPG) 

Video S2 Swimming behavior of Tg [smydl :m3ck) 
zebrafish at 13°C. The timeline represents: 00:00-00:26: After 
the fish moved into cold water bath and the beginning of cold 
treatment; 00:27-01:00: after cold treated for 1 h; 01:01-01:31: 
after cold treated for 2 h. 
(MPG) 
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